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Mild, unavoidable deviations from circular-symmetry of instrumental beams in current Cosmic 
Microwave Background (CMB) experiments now pose a significant challenge to deriving high preci- 
sion inferences from the high sensitivity and resolution of CMB measurements. It is important to be 
able to measure and characterize this subtle effect since it has bearing all subsequent inferences, in- 
cluding the cosmological parameters. We present analytic results, verified by numerical simulations, 
that CMB maps of cosmological signal that respect underlying statistical isotropy (SI) symmetry, 
measured with an instrument that has mildly non-circular (NC) beam would, nevertheless, exhibit 
SI violation. Further, we show that appropriate observable measures constructed within the Bipolar 
spherical harmonic (BipoSH) representation of SI violation capture subtle NC-beam effects coupled 
with the scan strategy of the instrument. 

Accompanying their latest 7- year data release, the WMAP team published very high significance 
measurements of non-zero BipoSH spectra, and A^°_ 2 ii m the "W" and "V" band of the exper- 
iment. The BipoSH measurements at the two frequency channels show significant differences that 
point against an origin in the achromatic cosmological CMB anisotropy signal. We present a strong 
case that the BipoSH measurement are primarily explained by the quadrupolar (m = 2) component 
of the NC-beams, b(2, of the respective channels. The fact that subtle levels of non-circularity, e.g., 
WMAP beams |fej2 |/&eo < 0.01, lead to measurable BipoSH spectra points to the immense promise 
and potential of the BipoSH representation. We recover non-trivial qualitative features of the mea- 
sured BipoSH spectrum that may be difficult to mimic through other effects. The scan strategy 
of CMB experiments, typically, employ multiple visits with varying orientation at different sky po- 
sitions. Hence, the effective smoothing associated with the NC part of the beam is not expected 
to be represented by the raw NC-beam map. In general, the effective beam at any point will also 
have variations over the sky. The key result of this work is that using BipoSH measurements it is 
possible to estimate an equivalent single hit, 'parallel transported' effective NC-beam that matches 
both the angular power spectrum and the non-zero BipoSH measurements of the observed maps. 
Hence, BipoSH analysis provides a very simple and effective characterization of CMB maps made 
with NC-beam. Last, but not the least, we successfully characterize an important systematic effect 
that limits all attempts to probe the SI assumption implied in the cosmological principle at high 
precision through CMB measurements beyond the angular power spectrum. 



I. INTRODUCTION 

Cosmic Microwave Background (CMB) experiments with high sensitivity, fine angular resolution covering almost 
the full sky have reached a point where even mild, unavoidable, deviations from circularity of the instrumental beam 
around its pointing direction poses a challenging systematic effect that must be overcome within finite computational 
resources to fully realize potential of these measurements for precision cosmology. 

The effect of non-circular (NC) beam on angular power spectrum of CMB has been studied in literature and the 
non-trivial impact on high precision cosmological inferences has been appreciated but not satisfactorily resolved, 
particularly, within available computational resources. A significant body of literature attempting to deal with NC 
beam effect on the angular power spectrum exists, E.g., [1-8]. Further, the beam response function can acquire an 
inherent time dependence which complicates and hinders the deconvolution enormously. Consequently, beam imper- 
fections, coupled with the scan strategy, lead to very complex modification of the signal demanding high computational 
resources to assess the final effect on the estimation of angular power spectrum and the cosmological parameters. 

It is important to be able to capture the NC-beam effect in observed CMB maps. We provide a computationally 
efficient characterization of the NC-beam effect by focusing on its effect beyond the angular power spectrum. Beyond 
the angular power spectrum, the improvement in quality of CMB measurements now allow observational confrontation 



'Electronic address: nidhijoshi@ctp-jamia.res.in 

^Electronic address: santanud@iucaa.ernet.in aditya@iucaa.ernet.in sanjit@iucaa.ernet.in tarun@iucaa.ernet.in 



2 



of simplifying, theoretical, assumptions implicit in current cosmology. Given CMB fluctuations are a realization of 
correlated, Gaussian random field, the two point correlation function uniquely describes the statistics of the fluctua- 
tions. SI implies rotational invariance of the two point correlation function. Breakdown of SI can be parametrized by 
the expansion coefficients of the two point correlation function in Bipolar spherical harmonic (BipoSH) basis known 
as BipoSH coefficients. Cosmological CMB temperature fluctuations are generally assumed to be a realization of SI, 
Gaussian, correlated, random field on the sphere, consequently, the angular power spectrum has been the primary 
observational target of most CMB experiments. However, current and upcoming CMB experiments also hold the 
promise to observationally constrain the underlying, often implicit, SI assumption (closely linked to the so called, 
'cosmological principle') employed in cosmology. 

Statistical isotropy (SI) assumption underlying modern cosmology has been under intense scrutiny with hints of 
various 'anomalies' persisting in successive years of WMAP data [9-13]. Violation of SI can be generated both, from 
theoretically motivated possibilities such as, cosmic topology, anisotropic cosmologies etc. as well as from observational 
artifacts, such as, beam non circularity, anisotropic noise, foreground residuals, etc. [14-23]. Searching for the most 
plausible candidate of SI violation is a non trivial endeavor which strongly depends on observational and theoretical 
hints to narrow down the the possibilities. While NC-beam effects do pose as a serious systematic contaminant for SI 
measurements, we demonstrate that the BipoSH representation measures of SI violation are also extremely effective 
in capturing, characterizing, and possibly, isolating subtle NC-beam effect in CMB maps. 

The recent detection of a quadrupolar power asymmetry 'anomaly' in WMAP 7 year data (WMAP-7) by the 
WMAP team using BipoSH representation measures [24] has proved to be an intriguing observation, that is yet to be 
satisfactorily explained. The published detections were significant both in the V-Band and W-Band frequency maps. 
Due to mildly significant differences in the BipoSH signal at the two frequencies and the fact that the signal seems 
to exhibit azimuthal symmetry in the ecliptic coordinate, it was suspected not to be of cosmological origin. Effects 
of residual galactic foreground emission would, arguably, be expected to be associated symmetries in the Galactic 
coordinate. 

Observed CMB anisotropy (and polarization) on the sky are a convolution of the underlying cosmological CMB 
signal with the instrumental beam response function. The instrumental beams in most CMB experiments are nearly 
circularly (azimuthal) symmetric, but mild NC deviations do arise due unavoidable limitations in instrumental design, 
function and fabrication, e.g., the primary lobe of the beam exhibits non circularity due to the off-axis position of 
detectors on the focal plane, diffraction around the edges of instrument leads to side lobes of the beam, due to finite 
response time of detectors the scan may not correspond to the instrument rotating around its beam axis leading to 
the effective beam response at any pointing direction being sensitive to the scan strategy, etc.. Regardless of the 
specific origin of non circularity in the beam response function, the key common point is the potentially measurable SI 
violation generated in the observed CMB maps. It has been argued whether the non-circular beam effect in WMAP 
provides a plausible explanation of the observed non-zero BipoSH measurements [25, 26]. This paper presents the 
first results from a research program to directly assess the nature and amplitude of SI violation that can arise due to 
NC-beam response within the BipoSH representation and its implication for the WMAP-7 BipoSH 'anomaly'. 

In the absence of symmetries, NC-Beam functions can be most generally expanded in BipoSH basis and coefficients 
of the expansion, Bf^ , are referred to as beam-BipoSH. An ideal, circular symmetric beam ensures vanishing beam- 
BipoSH for all L > 0. Importantly, the beam-BipoSH at L > not only capture NC-beam shape but also include the 
additional modulation arising from the specific scan strategy that sets the orientation of the NC-beam at any pixel. 
We show that every non-zero beam-BipoSH coefficient, iJj ^ , would generate a corresponding non-zero CMB BipoSH 
coefficient, AfYi in the observed CMB map. 

The analytic formulation presented is valid for any arbitrary beam shape but progress to explicit expressions is 
possible within an idealized 'parallel-transport' (PT) scan strategy approximation [1]. Further, mild deviations from 
circular symmetry, that permit a perturbative approach, and residual, discrete even-fold azimuthal, symmetry in the 
NC beam imply that the analytic results at the leading order quadrupolar-non-circularity (m = 2) are sufficient for 
the present observations. At the leading order of our approximations, we provide simple explicit analytic expressions 
only for BipoSH coefficients due to mildly NC-beams that retain discrete even-fold azimuthal and reflection symmetry 
relevant to explaining the WMAP-7 BipoSH results. Assuming reflection-symmetric beam function restricts the set to 
even parity BipoSH coefficients [45] . The published WMAP BipoSH results claim that in ecliptic coordinates BipoSH 
at (L = 2, M = 0) are the non-zero BipoSH measured. Our general analytic formulation then suggests that a PT- 
scan [1]) approximation should be fairly good simplifying approximation for the WMAP scan for the dominant m = 2 
mode of NC-beams in the ecliptic coordinates. We support this assertion with maps quantifying the departures from 
this approximation for a realistic WMAP scan. Further, the multiple hits of the beam with varying orientation at any 
pixel would tend to reduce (average out) the level of non-circularity significantly. In particular, it can be estimated 
in the specific case of WMAP scan, the m = 2 mode, bi2 , in the raw NC-beam can be expected to reduce by about a 
factor of ~ 0.5, and higher m modes maps are expected to suffer even progressively greater reduction. However, we 
also show analytically that scaling will have a dependence on the multipole, I, that is, in principle, computable given 
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detailed scan information including multiple visits and corresponding orientation at all the sky pixels. 

Assuming that the entire BipoSH signal measured by WMAP arises due to NC-beam effect, we determine spectrum 
of quadrupolar (m = 2) mode, bf^ , of an effective beam (incorporating multiple visits by the beam to sky pixels) that 
fits the measured BipoSH spectra. We fit for a constant rescaling of the raw beam = oibn . Although, the best- fit 
a ~ 0.4 obtained matches our expectation, the reduced x 2 °f the fit is unsatisfactory, and would, at best, explain part 
of the BipoSH signal. We find that a linear I dependent correction to a constant scaling, 6£f = (a + j3-l)bi2 lead to very 
good fit to the measured WMAP-7 BipoSH spectra. We provide analytic understanding of the /-dependent scaling, 
//. The linear fit explored is a simple 'phenomcnological' fit to Detailed, computationally expensive numerical 
simulations incorporating the two beam differencing, real scan and map making may justify the linear correction to 
the constant scaling, or, lead to the conclusion that the entire BipoSH signal cannot be attributed to NC-beam leaving 
room for other uncorrected systematic effects, or even hint at cosmological SI violation. 

We verify all analytic results and estimate error-bars through extensive numerical simulations of SI CMB maps 
scanned by real space convolution with corresponding NC-beam maps. The numerical simulations presented here are 
limited to the PT-scan, single side, NC-beam studied analytically. As explained later, without simplifying assumptions, 
realistic simulations using the WMAP beam differencing and scan are more compute intensive primarily because we 
find that the results arc fairly sensitive to side lobe and other spread-out features (out to ~ 10 times the FWHM 
from the beam centre), present in the beam maps released by the WMAP team [27]. Hence, we defer the results from 
ongoing realistic simulations to a future publication. 

The paper is arranged as follows. Sec. II provides a brief primer to the BipoSH formalism to characterize SI 
violations to keep the paper self contained. In Sec. Ill we present a novel expansion of the beam response function in 
the BipoSH basis, referred to as beam-BipoSH coefficients in the article. In Sec. IV we derive expressions for the CMB 
BipoSH coefficients arising from the convolution of SI CMB anisotropy signal with general NC-bcams. We also provide 
the simpler explicit analytic expressions for beam-BipoSH coefficients for a mildly non-circular beam within the PT- 
scan approximation. As a test case, we study the Elliptical Gaussian NC-beam model that has minimal parameters 
with clear interpretation. The beam spherical harmonic transforms (beam-SH) have analytic expressions and permit 
simple, well controlled perturbativc treatment of the NC-cffects. We also use this case to test and qualify our numerical 
simulations against analytic results. In Sec. V, we get closer to addressing the WMAP-7 BipoSH measurements by 
first computing the BipoSH spectra based on the beam-SH of single-side, raw beam maps (Sec. VA). In Sec. VB 
we fit for parametrized scaling functions, /; of the leading order raw beam-SH that would reproduce the WMAP-7 
BipoSH spectra measurements as arising entirely from uncorrected NC beam effect. Sec. VI has discussions and 
presents the conclusions of this paper. Detailed steps of all analytical calculations are provided for completeness in 
Appendix A and B. Appendix C provides an analytic derivation of an effective beam arising from multiple hits at 
varying orientations in terms of a rescaled beam-SH the raw beam map. 



II. PRIMER: BIPOLAR SPHERICAL HARMONIC (BipoSH) REPRESENTATION 

The CMB anisotropy signal is a random field on the surface of a sphere and can, hence, be expanded in the spherical 
harmonic (SH) basis, 

AT(n)=Y / ai m Y lm (h), (1) 

Vm 

and equivalently represented in the set of random variables et/ m . Statistical properties of a random field are charac- 
terized by its n-point correlation function. Statistical isotropy implies expectation values of all iV-point correlation 
functions are invariant under the rotations of the sky. The statistical distribution of fluctuations in temperature field 
are observationally consistent with Gaussian statistics. Hence, the two-point function C(h,h') = (AT(n)AT(n')) is 
expected to completely encode all information pertaining to the CMB anisotropy, AT(n). For CMB sky respecting 
statistically isotropy (SI), the correlation function only depends on the angle between the two directions h\ and iii 
and the two point function can then be expanded in Legendre polynomials, as 

C(n u h 2 ) = C(ni ■ n a ) = £ (2 * + ^ Cfci^ni ■ n 2 ), (2) 

l 

where the coefficients of expansion C\ define the angular power spectrum. The random spherical harmonic moments, 
ai m are then statistically independent leading to a diagonal covariancc matrix in the SH basis, 



{ai m a*i m ,) = CiSwSmm' 



(3) 
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which also implies that (the m-independent), C; encodes all the information in a Gaussian, SI field on the full sky 
(complete sphere, S 2 ). In general, SI violation leads to non-zero off-diagonal elements (for specific subset of SI 
violations, the SH covariance matrix remains diagonal, but m dependent). The two point correlation function, most 
generally, then depends on both the directions hi and h 2 and not just on the angle between them. The most convenient 
expansion basis for C(hi,h 2 ) ^ C{h\ ■ h%) is the Bipolar Spherical Harmonic(BipoSH) basis[28-33]. 

C(h 1 ,h 2 )= ^{Y h {h 1 )®Y h {h 2 )}LM, (4) 

h,l 2 ,L,M 

where A\™ are BipoSH coefficients and the bipolar spherical harmonic (BiPoSH) functions, 

{Y h (nx) ®Y l2 (n 2 )}LM = ^i^^i™!^) Y hril2 (h 2 ), (5) 

m\rri2 

are irreducible tensor product of two spherical harmonics spaces that form an orthonormal basis on S 2 x S 2 and, 
here, Cf^^i m2 are Clebsch-Gordon coefficients. The indexes of these coefficients satisfy the triangularity conditions 
\l\ — h\ < L < 1% + 12 and m\ + m 2 = M. The transformation properties of BipoSH under rotations are similar to 
spherical harmonics [34]. 

The BipoSH coefficients can be shown to be linear combinations of off-diagonal elements of the harmonic space 
covariance matrix, 

A-hh = ^ ( a hn l ia-l 2m2 )(^ 1 ) m2C lim 1 l 2 -rn 2 i ( 6 ) 
mim 2 

where a; m 's are the spherical harmonic coefficients of the CMB maps [28]. The SI part given in Eq. (3), corresponds 
to (L = 0,M = 0) , since , C^mi'-m' ^ ^ll'^mm'- Non-zero BipoSH coefficients with L > that capture SI violation also 
neatly encode the residual symmetries of the two point correlation function [3- r >] . The correlation function is always 
real and symmetric under the exchange of hi and h 2 , reflected in the following symmetry properties for BipoSH 
coefficients, 

{(_l\h+h-L+M £*L-M 
\ I hl 2 
(7) 
{ _ 1)h+ lz-L A LM 

As a consequence, BipoSH coefficients are always real if M = and l\ + l 2 + L is even and always imaginary if 
M = and l\ + l 2 + L is odd. If the observed CMB sky is statistically isotropic, it can be shown that all the BipoSH 
coefficients vanish except the coefficients of the form A^j which can be expressed in terms of the CMB angular power 
spectrum [28], 

Aff 2 = (-l) h C h l[6 hl2 6 L odMo, (8) 



where, for brevity, the standard notation Yi a b c = V + l)(2b + l)...(2c + 1) is introduced and employed in the 
rest of this paper [34]. 

Recently it has also been realized that BipoSH space representation separates into two distinct classes, 

• Even parity BipoSH coefficients Af^^ where l\ + l 2 + L is even obey Af^ (+) * = (—l) M Af i f M+ , and 

• Odd parity BipoSH coefficients Afjf - ' with h+l 2 + L odd, obey Afjf M * = (-1) M+1 Af^^ , 

where the 'parity' terminology here is based on the analogy with parity transformation properties of ordinary spherical 
harmonic moments. 

It is often easy to anticipate, or determine, the parity property of a possible cause of SI violation. This distinction 
then provides very valuable clues to the origin of SI violations e.g., weak lensing due to scalar (even parity) and 
tensor (odd parity) perturbations [36], anisotropic primordial power spectrum (even) [18], temperature modulation 
(even) [37], primordial homogeneous magnetic fields (even) [17, 38]. Importantly, in the context of origin of WMAP 
BipoSH measurements from NC-beam effect, the absence of significant odd parity BipoSH would justify the simple 
parallel-transport (PT) scan approximation adopted in our analytic approach. 
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III. Beam - BipoSH: NON-CIRCULAR BEAMS IN BipoSH REPRESENTATION 

It is very important to characterize the non-circularity of the beam through appropriate measures. We show that 
NC-beam is clearly represented in BipoSH space. Beam response function characterizes the angular dependence of 
the sensitivity of the apparatus around the pointing direction h\ . We expand the beam function around the pointing 
direction h\ in Spherical Harmonic (SH) basis as, 

B(hx,n 2 ) = ^h m (nx)Yi m (n 2 ). (9) 

l m 

The beam-SH coefficients when the beam is pointed at ni, are given by 

bi m .{hi) = J dil fl2 B(n l7 h 2 )Y hn (h 2 ). (10) 

The SH transform of beam at arbitrary pointing direction, n = (9, <f>) is given by beam-SH, bim 1 (z) - the SH transform 
of the beam pointed along z, the North pole [1], 

km(n) = J2 b lm>{z)D l mm ,{<j), 9, p(fl)), (11) 
m' 

where Wigner D-functions D l , (a, (3, 7) are the matrix elements of the rotation operator (0 < a < 2ir, < /3 < 
7r, < 7 < 27r) and a, /3,7 are the Euler angles that rotates the z-axis to the pointing direction h = (9, </>) and the 
angle p(n) specifies the orientation of the NC-beam with respect to the local Cartesian (x = 4>,y = 9) aligned with 
spherical coordinates [1]. Such a rotation can be realized by fixing a coordinate system and performing anti-clockwise 
rotations, first rotating about the z-axis by an angle a = </>, then rotating about new y-axis by an angle (3 = 9, and 
finally about the new z-axis by 7 = p(n). Rotation is such that the beam pointing direction, n and p(n) is the relative 
orientation of the beam around its axis set by the the scan pattern of the experiment . 

In this work we introduce a more general NC-beam characterization in the BipoSH representation. Since, a general 
NC-Beam function depends on two vector directions, it can be expanded in the BipoSH basis (see Sec. II) , 

B(n 1 ,n 2 )= £ Bffi Yl Ch^rnJhrn^YlzrnAn^ (12) 
Ix^LM m 1 m2 

where the coefficients of expansion B^ are referred to as beam-BipoSH coefficients . As will be clear in the sections 
ahead, beam-BipoSH coefficients have the advantage of capturing the additional effect of scan strategy together with 
non-circularity of the beam. 

The beam-BipoSH coefficient can be readily related to beam function in spherical harmonic basis , 

B h% = E C M 2 m 2 I dn A b l2m2 (n)Y* mi (n) . (13) 

mini2 

A circularly symmetric beam function around the pointing direction can be expanded in Legendre polynomials, 
B(hi,n 2 ) = B(fi\ ■ h 2 ) = (47r) -1 ^ ; (2Z + l)BxPi{h\ ■ n 2 ). Inverse transforming Eq.(12) and using orthogonality of 
BipoSH [34], we obtain beam-BipoSH coefficients for circularly symmetric beam function, 

B?£ = B h S hl2 (-l) h 1[[S L0 S M0 , (14) 

h 

where Bi is the commonly used Legendre transform of the beam function in the circularized beam approximation. 



A. General expression for an arbitrary scan strategy 



Using Eq. (11) and (13), it turns out that for any arbitrary scanning strategy, the beam-BipoSH can be expressed 
in terms of the beam-SH as, 

B hh= E C&ar^E^'C 8 ) r / " D l ^ m ,(<t>,6,p(6,<t>))Y£ mi (6,cj>)sin6d6dct>. (15) 

„ J9=0J<t>=0 
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To separate the azimuthal (0) and polar (9) dependencies, it is convenient to expressed Wigner-D functions in terms 
of Wigner-d through following relation, 

D l mm , (0, 9, p) = e-*»* 4m' (») c "" n ' P - (16) 

Eq. (15) is the most general expression of beam-BipoSH coefheicnts for any given NC-beam specified through bi m (z) 
and scan pattern, defined by p(0,<p), in any spherical polar coordinate system. The beam-BipoSH coefficients can 
be evaluated numerically from the beam maps of the NC-beam. When the beam is circularly symmetric, bi m (z) = 

Bi J ,j to0 ; anc ] it i s straightforward to establish that Eq.(15) reduces to Eq.(14). Beam-BipoSH depend not only 

on NC-beam harmonics but also on the scan-strategy that defines p(n). We motivate a particular idealized scan 
pattern where p(h) is constant in a some coordinate system, that not only allows a completely analytic treatment 
but is also reasonably well-justified in the context of the BipoSH signature measured in WMAP-7. Analytic progress 
is far less tedious when the beam has mild deviations from circularity allowing a perturbative approach that retains 
only the leading order terms. We validate numerically that these approximations used are indeed fairly adequate for 
the explanation of the WMAP BipoSH measurements. 



B. 'Parallel-transport' scan approximation 

The general beam-BipoSH in Eq. (15) can be tackled analytically if when the scan pattern is such that p(h) is a 
constant. We refer to this case as 'parallel-transport' (PT) scan following [lj. The PT-scan approximation implies 
that the orientation of the beam relative to the local longitude is constant at any point on the sky. A beam response 
function with reflection symmetry will generate only real, non-vanishing beam-BipoSH coefficients. The leading order 
m = 2 dominates for a mildly NC-beam. For real coefficients only the cos(2p) part is important (see Eq.(15)). Fig. 1 
shows the map for cos(2p) for a full year scan. For WMAP scan pattern, the constant p approximation for m = 2 
holds good over the significant (blue) band around the equator in an ecliptic coordinate system. Therefore, for a 
fair fraction of the sky part of the map our constant p approximation for m = 2 mode of non-circularity is a fair 
assumption. The compute-intensive numerical comparison of results from this approximation with the real WMAP 
scan is underway and will be reported in the near future. 

For constant p(n) = p, the orthogonality relation, 

/ exp(-«(mi + ma)0) defy = 2ir 6 mi ,-m 2 , (17) 
J<t>=o 

implies that the integral over <p i n Eq.(15), will separate from the integral over 9 and would restrict the non-zero 
beam-BipoSH to M = 0, 

B™ = 27^0^^^ £*W (*) expHm'p) £ ff£° maZam2 f d l ^ m ,(9)d l l m20 (9) sin Odd. (18) 

m' m 2 J 8—0 

We define 

Ittm' = (-^V^F / ^ 2 m'W4 2 o(0)sinM?. (19) 

Here we have used the symmetry property of Wigner-d functions, 4m' = (— l) m_m d_ m _ m ,. Beam-BipoSH coeffi- 
cients in the constant p(n) approximation take the following form, 

Bffi = 2nS M0 k 2m ,(z) eM-im'p) £ C™ m2l2m2 I^ m , . (20) 
m' »2 

For a mildly non circular beam, the summation can be truncated at a low value of to'. Further, NC-beams that retains 
discrete even- fold azimuthal symmetry will have bi m (z) = for odd values of m. In BipoSH space, the consequence 
of discrete even-fold azimuthal and reflection symmetric NC-beam translates to restricting non-zero beam-BipoSH to 
AI = and l\ + 1% = even. Imposing reflection symmetry, beam-BipoSH are restricted to even values of multipole L. 
We derive explicit expressions for beam-BipoSH coefficients up to leading order (w! = 2) term . 



FIG. 1: Top: Map of p(h) = [A(n)] _1 J2i Pii™) defining the orientation of the beam averaged over multiple hits, 
i = 1 to N(n) for numerically simulated one year of WMAP scan. At most pixels, p(n) averages to a value close to 

zero. However, it also shows a very intricate fine patterns of small regions with marked deviation from zero. 
Bottom: Map of cos(2,o) in ecliptic coordinates for full year WMAP scan provides a rough guide to the validity of 
the PT-scan approximation for the dominant m = 2 non-circularity for the WMAP scan strategy. Note that over a 
wide band of latitudes (blue region), the PT-scan approximation is roughly valid. Strong deviation are seen at and 
beyond the two latitude-symmetric thin bands (red) at moderately high latitudes in North and South hemispheres. 

To make analytical progress, we need to evaluate I^m' f° r m ' = — ^, 0, 2. It is important to note that, m' = 
corresponds to the usual circular symmetric beam-BipoSH (L = 0) coefficient. The non-circular m! = ±2 will give rise 
to non-trivial (quadrupolar, L = 2) beam-BipoSH. For circular symmetric (to' = 0) beam-BipoSH (refer AppcndixA), 




(21) 



,LM 

'hh 



{-l) l2 bi 20 (z)V^7rSi 1 i 2 S L0 5MoS m 'o- 



(22) 



For m' = ±2, the integrals are evaluated separately for the = and mi ^ parts of the summation. In the former 
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case when mi = and ml — ±2, the integral 



rhh 

i 0,±2 



(2/i + 1) 

47T 



. (h-2)\ 



C2-2)! r 4ig 
(Z 2 +2)! L(2i 2 + 1) 



2i 2 (i 2 + l) - 
(2i 2 + l) - 



if (/! + ? 2 = Odd) 

if (Zi > Ja) 
if (li < i a ) 

if (h = h) ■ 

hh 



(23) 



For m 2 ^ 0, dj^ ± 2 (Q) i s recursively expanded in terms of d£ o (0) to evaluate 1^ ±2 (refer Appendix A). 

This work is motivated by the highly significant measurements of BipoSH published by the WMAP team [24]. 
Interestingly, the publication claimed that for the L = 2 BipoSH spectra measured are significantly non-zero only 
for M = in the ecliptic coordinates. This has been independently confirmed by our BipoSH measurements on 
WMAP-7, as well. As shown in following Sec. IV, this implies that only beam-BiPoSH with M = are relevant for 
understanding the WMAP BipoSH results. Consequently, it suggests that if NC-beam is responsible for non-zero 
BipoSH measurements, the WMAP scan pattern is such that, in ecliptic coordinate, p(h) w constant must be a good 
approximation for the dominant m = 2 mode of the NC beams of WMAP. This is being quantitatively assessed through 
more compute intensive, numerical simulations with realistic WMAP scan patterns (i.e., without invoking PT-scan 
approximation and incorporating map- making from actual two side difference measurement). We take advantage of 
the simplicity afforded in the PT-scan approximation to proceed to explicit analytic expressions that readily provide 
broader as well deeper insight than compute expensive numerical simulations. 

This also implies that only even-parity beam-BipoSH would be relevant for our comparison to WMAP BipoSH 
results where the estimator used is by definition is strictly restricted even-parity BipoSH. Non circular beam with 
additional reflection symmetry restrict beam-BipoSH to even-parity, 



B, 



LM _ niM (+) 



hh 



= B 



hh 



= 2n5 MO (b h2 (z)e- i2 " + b] 



(S)e i2 ') 



/~iLQ jhh 

°;iO/ 2 o J o,2 



E r L0 j 
^h-m 2 h™2 • 



hh 
m 2 ,2 



m 2 #0 



(24) 



Refer to Appendix A for details. Note that a constant p can be absorbed as phase factor in the redefinition of the 
complex quantity 6; m (i) essentially resetting the orientation of the beam when pointed at North pole. Hence, in the 
PT-scan approximation, 
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(25) 



We emphasize that the above expression for beam-BipoSH coefficient hold for the PT-scan approximation (with 
constant p(nj) for a NC-beam that has reflection symmetry. Although we have restricted explicit analytic results 
presented in the text to reflection symmetric beam functions, in general, odd parity beam BipoSH will be non- 
vanishing in absence of the above mentioned symmetries. Appendix A provides expressions for odd-Parity beam- 
BipoSH, B££f ( ~\ that can be used as a measure of breakdown of reflection symmetry in NC beam[4G]. Note that 
the BipoSH estimator [37], that differ by a factor from original definition of Hajian & Souradeep [28, 33], used by the 
WMAP team cannot be extended to odd-parity BipoSH, However, it is possible to devise BipoSH estimators that can 
measure odd-parity BipoSH spectra while matching that employed by WMAP for even-parity BipoSH spectra [39]. 



IV. NON-CIRCULAR BEAM IMPRINT ON BipoSH OF CMB MAPS 

The observed CMB map is a convolution of the signal with the instrumental beam. Instrument beam response 
functions of all experiments have deviations from circular symmetry at some level. However, circular symmetric beam 
response function around the pointing direction is often assumed to simplify the analysis of the the beam effect. This 
key assumption does affect many stages of CMB data analysis and should be measured and characterized. 

In this section we show that CMB maps made with non-circular beam exhibit statistical isotropy (SI) violation in 
an otherwise SI cosmological CMB signal. BipoSH coefficients have proved to be robust, model independent, measure 
of SI violation. Further, BipoSH representation provides clear insight into the nature and extent of the residual 
symmetry and parity [35, 36] and, hence, helps pin down the most plausible origin of SI violation. 
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The measured CMB temperature fluctuations map, AT(n) is the convolution of true underlying CMB temperature 
fluctuations AT(n) with the instrument beam, 



AT(ni) = / <m n2 B(n u h 2 )AT(n 2 ). (26) 

where the beam response function B(ni, h 2 ) gives the sensitivity of an experiment around the pointing direction. The 
observed two point correlation function will be, 

C{h u h 2 ) = (AT(ni)AT(n 2 )) = / dCt n [ dQ„'C(h' ,h)B(ni,n')B(n 2 ,n), (27) 



where C(n',n) = (AT(h')AT(n)} is the underlying correlation function relevant for cosmology. It is evident from 
Eq.(27), that SI violation can occur either due rotational invariance breakdown of the true underlying temperature 
correlation function, C(n 1 ,n 2 ) i C(hi ■ n 2 ), or, due to the breakdown of circularity in beam response function 
B(n ll h 2 ) it B(h 1 ■ n 2 ). 

If the underlying CMB signal respects statistical isotropy (SI) symmetry as widely assumed in cosmology, then 

9/ i 1 

C(n lf n 2 ) =J2^—Q W^huhs), (28) 
i 

where Wi{h\, n 2 ) is the elementary window function [I] that accounts for effect of the finite resolution of beam function, 
given as 



Wi{hx,n 2 ) = J dft fi J dSh'B(fii t b)B(th,n')Pi(h-ti). (29) 
Inverse transform of Eq.(4), yields the CMB BipoSH coefficients, Af 1 ^ in terms of beam-BipoSH coefhcients as 

Affi = [ dn ni f dQ^CifhMYiM^Y^na)}^ (30) 



Hence, when the underlying CMB signal is statistically isotropic and beams are circular, BipoSH coefficient are 
expected to vanish for L > 0, refer Eq.(8), 

Aff 2 = (~l) h l[C h BlS lll2 6 LO SMo. (31) 

h 

In Sec. Ill, we show that non-circularity of the beam response function is captured by beam-BipoSH coefficients. 
The BipoSH coefficients of CMB map for a reflection symmetric beam function can then be derived in terms of 
beam-BipoSH coefficients as (see Appendix A), 



A hh ~ l^y~ l ) Ll 2^ B hl B hl LL C L 1 M 1 L 2 M 2 \ L T l 
i L 1 M 1 L 2 M 2 L\L" 



(32) 



In Sec. Ill, we also argued that it is convenient, as well as observationally motivated, to carry out the analytic analysis 
in a coordinate system where parallel-transport (PT) scan approximation holds and, consequently, the non-zero beam- 
BipoSH are restricted to M = 0. Eq. (32) then dictates that the corresponding BipoSH coefficients of the CMB maps 
are also restricted to M = and are given by 



ki = z^m e n ^v 2 o{l l L \ t 2 2 } • (33) 

Since there is a triangularity condition (\L\ — L 2 \ < L < L\ + L 2 ), therefore the most dominant term in the above 
summation are {Li = L,L 2 — 0} and {Li = 0, L 2 = L}. 

The BipoSH estimator used by the WMAP team [24, 37], that differs from the original BipoSH definition in Hajian 
& Souradeep [28] by a factor of IIl/CILi!, ^hoho) are restricted to only even-parity BipoSH. We insert this factor 
in the BipoSH expression in order to allow comparison to the published WMAP measurements more transparent [47] 
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(34) 



SI violation signals in WMAP-7 were measured in two BipoSH spectra, Af® and Af® 2[ , we provide explicit leading 
order expressions for these coefficients arising from the NC-beam as, 



(Ui) 3 Cf % 



(35) 
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nv 



n, 



(36) 



Note, that BipoSH expression in Eqs. (35) and (36) are provided in the scaled form that matches the BipoSH estimator 
employed by the WMAP team. 

A. BipoSH from test Elliptical Gaussian beam model 

Elliptic-Gaussian (EG) functions provide a simple model NC-beam extension to the oft-used circular-symmetric 
Gaussian beam function. The non-circularity is clearly parametrized by the eccentricity of the elliptical iso-contour 
lines. An EG-bcam function pointed along z axis can be expressed in spherical polar coordinates, as 



B(z,n) 



1 



27r<7l(72 



exp 



(37) 



where the azimuth angle dependent beam- width <j{4>\) = [af/(l + e sin 2 t/fi)] 1 / 2 is given by Gaussian widths o\ 
and (72 along the semi-major and semi-minor axes. The non-circularity parameter e = (cr 2 /u 2 — 1) is related to 
eccentricity e = y/l — a 2 ja\ of the elliptical iso-contours. As expected, the Elliptical-Gaussian beam reduces to 
circular Gaussian beam for zero eccentricity (e = 0). Higher the value of eccentricity, stronger the deviation from 
circularity. An analytical expression for the bcam-SH of EG-beam is available 



hm{z) 



(2l + l)(l+m)\ (] ly 
47r(i-m)! V ^ 2 > 



L m/2 



exp [— {I 



2> 2 



-)] (m = even) 



(m = odd) 



(38) 



where I v {x) is the modified Besscl function [1]. The discrete even- fold azimuthal symmetry and reflection symmetry 
dictates bi m = for odd to. The reality condition of beam, b* m = bi m for even to, then implies 6;_ TO = b\ m . 

For Elliptical-Gaussian beams, the ratio bi m /bio dies down rapidly with |to|. Hence, for mild eccentricities of the 
Elliptical beam, it is sufficient to consider only the lowest to = 2 mode. In Fig. 2, we plot spherical harmonic transform 
of an Elliptical-Gaussian beam with Qfwhm = 13.579', which is close to Qfwhm of W band of WMAP (with ellipticity 
e = 0.4 that best matches WMAP-7 W-band BipoSH spectra amplitude). Similarly, we also consider a EG-beam 
with Qfwhm = 17.7036' close to the V-band beam. (The best fit to V-band BipoSH spectra corresponds ellipticity 
e = 0.46). For analysis purpose, it is a good enough approximation to restrict to m = 0, ±2 modes. We estimate 
beam-BipoSH coefficient in Eq.(25) for the case of "non-rotating" beams, p = by using the closed analytical form of 
bimS. Finally, using Eq.(35) and Eq.(36) we obtain the CMB BipoSH spectra Af® and A 20 _ 2l that were measured to 
be non-zero in WMAP-7. We verify the analytical results and compute the error bars on the BipoSH coefficients using 
statistically isotropic Gaussian simulations convolved with EG-beam functions. In Fig. 3, the analytic calculations of 
CMB-BipoSH for two cases of EG-beam are shown overlaid on estimates from average measurements with error-bars 
from 100 numerical simulations of noise free SI maps convolved with corresponding NC-beams. The results with test 
elliptical-Gaussian beams serve to crosscheck and validate our analytical expression derived and a check the numerical 
simulation of CMB maps convolved (in real space) with an NC-beam. 

It is clear that NC-beams generate detectable levels of non-zero CMB BipoSH. The CMB-BipoSH spectra Af® 
and Af® 2l are both proportional to the the leading order non-circular correction, biob^ and the underlying SI CMB 
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FIG. 2: Harmonic transforms of an elliptical Gaussian beam with Qfwhm = 13.579' and ellipticity e = 0.4. The 
circular symmetric component of the NC EG-beam is given by is to = and to = 2 captures the leading order 

EG-beam effect. 



power spectrum (see Eqs. (33) and (24)). It is instructive to compare the EG-model beam predictions to the WMAP 
BipoSH measurements. Using Elliptical-Gaussian beams with FWHM values close to the WMAP W and V band 
reproduces detectable peak structure in the BipoSH spectra. The amplitudes can be matched by choosing appropriate 
values of effective eccentricity for the V (e = 0.46) and W bands (e = 0.4). Although, the eccentricity of the WMAP 
beams can be estimated for raw beam maps [2], it is expected that multiple visits by the beam at any pixel with 
varying orientation in WMAP scan would tend to reduce the effective eccentricity to smaller values. The eccentricities 
that reproduce the correct amplitudes for the CMB-BipoSH roughly correspond to the reduction also indicated from 
comparing the corrections to the angular power spectrum [40]). The fact that both the Af® and Af^_ 2l spectra, are 
roughly matched in amplitude and qualitative peak structure with the correct relative sign, for identical ellipticity 
parameters, then suggests that the WMAP-7 BipoSH measurements may have actually captured subtle uncorrected 
NC-beam effects. 



V. BIPOSH SIGNATURES OF NON-CIRCULARITY IN WMAP BEAMS 

While CMB BipoSH from the test elliptical-Gaussian NC-beam motivates a more careful estimation of the NC- 
bcam effect in WMAP-7 maps, these fail to reproduce a key qualitative feature of WMAP BipoSH measurements; 
viz. the change in sign of the V-band Af® at large multipoles. This, however, is not surprising. It is widely- 
known that the WMAP beams deviate from a Gaussian profile and that an Elliptical Gaussian beam is not a good 
approximation [2, 41, 42]. In section VA, we first use the SH transform of the raw WMAP beam maps to compute 
the CMB-BipoSH under PT-scan approximation. Interestingly, we recover qualitative features in BipoSH spectra 
measurements that were not captured by the Elliptical Gaussian beams. As expected the amplitude of the CMB- 
BipoSH based raw beam are much higher in amplitude relative to the measurements. A rudimentary numerical 
superposition of raw beam transforms to mimic the effect of multiple hits with varying orientation at a pixel (using 
a WMAP scan), leads us to expect a constant scaling of the leading order bi2 by a factor a ~ 0.4. More a detailed 
analytics treatment presented in the Appendix C that allows for deviations from PT-scan approximation in averaging 
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FIG. 3: BipoSH spectra Af® and Af® 2l for an elliptical Gaussian beam with : (Top) Bfwhm = 13.579' and ellipticity 
e = 0.4; (Bottom) Ofwhm = 17.7036' and ellipticity e = 0.46 smooth (blue) curve are binned CMB-BipoSH 
computed using the analytic expressions. The red curve with corresponding error-bars are obtained from 100 
numerical simulated SI maps convolved with the same elliptical Gaussian beam. The underlying SI signal 
corresponds to the best fit A-CDM model from WMAP7 



over multiple hits, reveals multipole I dependent correction, /; to the constant scaling that depends on the detailed 
scan strategy of the CMB experiment over all multiple hits. 

The analytic expression for /; in Appendix C, for simplicity, pertains only to a single beam scan case. In particular, 
the approach cannot account for any complexities that may be introduced to the effective NC-beam that is used 
estimates the BipoSH spectra arising due to linear algebraic map-making procedure from two beam differencing. 
However, it motivates the exploration of phenomenological effective beam bf^ = Ylm' Flm'blm' expressed in terms 
of the beam-SH obtained from the raw beam maps. We proceed with the assumption that the WMAP BipoSH 
measurements arise from NC-beam effect and take the 'phcnomenologicar approach to determine b^ from BipoSH 
measurements. Limiting attention to the leading order contribution at m = 2 we parametrize the function /; = Fa. 
and fit to the measured BipoSH spectra. In section VB, we first fit for a constant value a that most closely reproduces 
the WMAP BipoSH measurements. While we find that a ~ 0.4 is the best-fit value, the goodness of fit is not very 
satisfactory. We find that a linear correction to the constant scaling, bf^ = (a + f3 ■ l)bi2 provides very good fits to 
the WMAP BipoSH measurements. Numerical simulations incorporating complexities of two beam differencing and 
map-making procedure and deviations from PT-scan approximation are needed to decide whether the parameters of 
the phenomenological scaling obtained are justified. 
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A. BipoSH coefficients from single-side WMAP raw beam map 

To mimic NC-beams closer to WMAP, we consider the A side raw beam maps of the V2 and Wl differencing 
Assembly (DA) of WMAP as representative of the V and W band beams, respectively. We compute the beam-SH 
coefficients for these numerically for use in semi-analytic estimate of the CMB BipoSH coefficients, using Eqs. (35) 
and (35). Fig. 4 gives plots of the circular, bio and leading order ba beam-SH coefficients of the W1A and V2A raw 
beam maps. Note that the 6/2 spectrum changes sign and takes negative values at large I - a key qualitative feature 
of the WMAP beams that cannot be captured in Elliptical-Gaussian models where bi2 does not change sign with /. 
The origin of this curious feature becomes apparent in the raw beam maps of the WlA and V2A channel shown in 
Fig. 5. The central part of the beam maps show an elliptical peak with marked non-trivial non-circular 'shoulder-like' 
features. However, more interesting are the right-hand panels, where we highlight a spread out annular distribution 
of regions with negative response. The integrated power in the negative beam response ~ 0.5 of the total power 
and has significant impact on the beam-SH. In particular, it leads to negative values of bi2 at I > 2cr h ~ 1 . We find 
that this feature is critical in recovering the qualitative feature in V-band Af® changing sign at large I seen in the 
WMAP measurements. Since such an unique correspondence between a beam-SH feature and the consequent CMB 
BipoSH is unlikely to be mimicked by other effects. We claim that this provides a strong hint that WMAP BipoSH 
measurements are linked to the NC-beam effect. 

We use PT-scan approximation in ecliptic coordinates given the fact that the WMAP BipoSH results indicate 
azimuthal symmetry of SI violation (M = 0) in ecliptic coordinates. The analytical results are verified using numerical 
simulations of the BipoSH coefficients Aj® and Aj® 2 i are shown in Fig. 6. As seen in Fig. 6, the amplitudes of the 
BipoSH coefficients do not match with WMAP7 year's observed detections. This is expected as the BipoSH coefficients 
are obtained from the A side beam of a single channel in W and V band without accounting for circularizing of the 
beam due to multiple hits with varying orientations. We collate and list some key similarity of these results to the 
detections of SI violation in WMAP-7 data: 

1 . As expected from analytic understanding, the numerical simulations also confirm that the most significant con- 
tribution is the M = mode in BipoSH spectra Af t and Af^ 2[ when the analysis is done in ecliptic coordinates, 
giving strength to our PT-scan approximation. In a coordinate system where PT-scan approximation is valid 
only M = mode should be significant. 

2. We notice the NC beam effect is larger in W band than in V band explaining the difference in detected SI 
violation signal at the two frequencies. 

3. We recover the change of sign of the BipoSH A£° measurements at large I in the V-band. 

4. The BipoSH coefficients from NC beam shows a prominent bump roughly around the first acoustic peak (I = 220) 
for both W and V band. This corresponds mainly to the scale picked by the underlying angular power spectrum 
C[. However, the precise peak location also depends on the peak in 6;o^(2 for each band and can account for 
differences in the peak location in the two bands shown in Figure 4. 

The exercise of computing BipoSH from single side raw beam maps under PT-scan approximation docs provide mo- 
tivation for a deeper analysis of the NC-beam effect in the WMAP-7 maps. Numerically simulated maps incorporating 
the two side WMAP beam differencing, scan and map-making details should settle the extent to which approxima- 
tions used here are valid and clearly settle the question of whether the entire WMAP-7 BipoSH signal arises as an 
effect of uncorrected NC-beam effects. Implementing numerical simulation with actual WMAP beam maps present 
serious computational challenge. Besides the main central peak having significant structure, the WMAP beams have 
a spread-out annular pattern of negative response that cannot be neglected. This implies that for a careful BipoSH 
analysis of SI violation, a convolution with much larger regions of the WMAP beam maps of W and V band, extending 
roughly to, ~ WOfwhm , is required making it computationally expensive. The beam maps of WlA and V2A channel 
are shown in Fig. 5. In particular, we highlight the spread-out negative part of the beam response that explains some 
unique qualitative features of the WMAP BipoSH measurements. The power in the negative part of the beam is more 
than 50 percent of its central positive peak and hence, non negligible. Work is in progress and results will be presented 
in the near future. 



B. Fitting BipoSH spectra measurements to an effective NC-beam 

In the analytic approach we show that the BipoSH coefficients scales as Af^ cx bi ± obi 2 2 for mild deviations from 
circular beams under PT-scan as given in Eq. (33). However, realistically, complex details of relative orientation of A 
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FIG. 4: Top: Beam spherical harmonic transforms, bio and bn of beam maps of A side of Wl and V2 DA. It is 
interesting to note that \ba\/biQ < 0.01 implying the BipoSH representation is sensitive to really really subtle levels 

of non-circularity in the beams. 
Bottom: The plot of the NC-beam leading order NC beam perturbation parameter 6;o&;2 vs. I . The solid lines 
correspond to beam-SH of the raw beam maps, the dashed lines show the b^b^ of the phenomenological best-fit 
effective beam-SH b^ = fib\% obtained by parametrized linear fits to /; using the the two measured BipoSH spectra, 
Af® and Af® 2l The plots also show that although BipoSH peak structure is largely set by the underlying angular 
power spectrum C'i of SI cosmological mode, small differences observed at the two different frequencies can arise 

because of the difference in the shapes of &zo&Z2- 



and B side of beams for each DA, the differencing scheme and map-making, the varying orientations for the multiple 
hits at each pixel in the actual scan strategy, etc., all could have bearing at finer levels on CMB BipoSH generated 
by the NC-beam. 

As shown in Appendix C, the averaging of the NC-beam at any pixel due to multiple-hits at different orientations 
can be expected to lead to simply a scaling of = fibii. This would correspondingly scale BipoSH spectra obtained 
from A side beam of Wl and V2 channel by // . 

Retaining PT scan approximation, the superposition of raw beam transforms leads us to expect a constant scaling 
/; = a. The value a ~ 0.4 most closely reproduces the WMAP BipoSH spectra. We confirm this to be consistent 
with numerical estimates from beam-SH superposition using pi(n) obtained from WMAP-likc multiple hits but with 
PT-scan approximation (see Appendix C. The numerical value of a depends on the unknown relative orientations of 
the A and B side beams (assumed to be identical) and matches the fit value of a for relative orientation ~ 140°. 

However, the goodness of fit for constant fi is not very satisfactory in explaining the entire WMAP-7 BipoSH 
spectra in terms of NC-beam systcmatics and suggests that /-dependence of /; may be important for a finer match. 
While, in principle, it is possible to determine /; for the WMAP scan by numerical superposition of beam-SH using 



15 




FIG. 5: Beam response function of A side of Wl differencing assembly (Top) and A side of V2 differencing assembly 
(Bottom) are shown. The Left-hand panels show 3D zoomed view of the central part of the beam shows a truncated 
view of the central beam peak to clearly highlight the elliptical contour. Marked non-circular shoulder features are 
seen in the central peak. The Right-hand panels cover the entire beam-map images to show the spread-out annular 
distribution of regions with negative response. (The red circles marks the central beam peak region). The power in 
negative response is ~ 0.5 of the positive power in the central peaks. The annular negative power distribution shows 
quadrupolar feature that modifies the beam-SH bi2 to take negative values at high / for the V-band (see Fig 6). The 
corresponding impact of this beam-SH is significant and crucial for understanding the intriguing qualitative 
zero-crossing feature in the V-band BipoSH spectra. It is apparent then that correctly accounting for WMAP 
NC-bcam effects numerically in our ongoing analysis, requires the convolution of almost the entire beam map region 
with the SI sky- map leading to enormous increase in computing costs (relative to using only the central peak). 

Eq. (Cll). Practically, it will be plagued by ambiguity of relative orientation of the A and B side NC-beams and 
many other finer details. Hence, we prefer to take a phenomenological approach and explore parametrized /; that 
fit the BipoSH measurements. Minimally, this provides an effective beam that would explain the observed WMAP-7 
BipoSH spectra under PT-scan approximation. We believe that this provides a computationally inexpensive approach 
to quantify the WMAP NC-beam systematic effect. 

The obvious step beyond constant scaling, would be to add a linear correction /; = al + /3 * I. As shown in Fig. 7) 
we find very good fits to the WMAP BipoSH measurements with this scaling. Interestingly enough, it seems possible 
to find values of a and j3 that simultaneously provide fairly good fits to both Aj® and Af® 2l BipoSH spectra. Note 
that the product biob^ still tends to zero for large / as seen in the bottom panel of Fig. 4. 
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FIG. 6: BipoSH spectra, Af® (Left) and Af® 2 i (Right) obtained for the raw beam maps A side of V2 channel ( Top) 

and A side of Wl channel (Bottom) of the WMAP experiment. Analytically evaluated BipoSH spectra (blue) 
overlaid on the average BipoSH spectra (red) with error bars obtained from 100 simulations of statistically isotropic 
CMB sky convolved Wl and V2 channel of WMAP. While the BipoSH spectra from raw beam-SH recover 
interesting qualitative features of the measure spectra, the amplitudes need to be corrected for the circularizing 
effect of averaging raw beam maps with multiple hits with varying orientation at any pixel. 



VI. DISCUSSIONS & CONCLUSION 



The observed CMB sky is a convolution of the cosmological signal with the instrumental beam response function 
of the experiment (assuming perfect removal of foreground emissions, etc.). The deconvolution of the beam effect 
from the signal is relatively straightforward for an ideal circular symmetric beam. Non-Circular (NC) deviations of 
the beam , however mild, are practically inevitable in all experiments, and affect the results obtained at the limits 
of the sensitivity and resolution of the recent experiments. The effect of WMAP NC-bcam on the angular power 
spectrum C; is not very pronounced and is significant only at very large multipoles. It is, hence, very interesting 
that CMB maps obtained with NC-beams disrupt the rotational invariance of the two point correlation function of 
the cosmological CMB signal, theoretically expected from the 'cosmological principle', leading to clearly measurable 
signatures of Statistical Isotropy (SI) violation. 

We show that SI violation measures in the Bipolar Spherical Harmonic (BipoSH) representation of the observed 
CMB maps, are well-suited to capture systematic NC-beam effect. It is important to note that even though the level 
of non circularity is WMAP beams is at ~ 1% percent, the BipoSH spectra generated by this effect are measurable. 
This points to immense promise and potential of the BipoSH representation also as a diagnostic tool for current and 
future CMB experiments. A key objective is also to assess whether the recent measurement of non-zero BipoSH 
spectra, Af® and Af°_ 2l in WMAP-7 in the V & W band maps [24] could arise from uncorrected NC-beam effects in 
WMAP-7 maps. 
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FIG. 7: Comparison of BipoSH spectra generated by effective NC-beams at the V and W bands to the published 
WMAP-7 BipoSH spectra measurements (approximately digitized from the Fig 16 of Ref. [24]). As described in 
Appendix C, an effective NC-beam can be constructed from from the raw beam-SH, bi m given by 
buL = X)m' Fim'bim'- The measured BipoSH spectra, AfP and Af® 2[ , are sourced at the leading order by the m = 2 
term. Hence we estimate the parametrized fits to the function /; = F\ 2 using the W and V band BipoSH 
measurements. We present the best fit results for a constant fi=a and a linear function fi = (a + j3 ■ I). The 
respective values of the reduced \ 2 of the fit and the best- fit parameters are provided in the figure. The constant, 
provides a fit that partially accounts for the BipoSH spectra detections in terms of WMAP NC-beam effect. 
However, remarkably an effective from the linear form of /; can satisfactorily explain the non-zero BipoSH 
spectra detected in WMAP-7 in terms of the corresponding V and W band effective NC-beams. 



In this paper, we introduce the novel and useful concept of expanding the NC-beam response function in the 
BipoSH basis to define beam-BipoSH coefficients. The beam-BipoSH not only incorporate the effect of the NC-beam 
(expressed in terms of the beam-SH - the spherical harmonic transform of the beam bi m (z) pointing at z), but also the 
additional effect of the scan strategy that determines the beam orientation angle, p{n) relative to the local longitude 
at any pixel. Neat and completely analytic results for CMB BipoSH generated by NC-beams can be obtained within 
the 'parallel-transport' (PT) scan approximation where the beam visits the pixels at constant p(n) = p. We argue 
that this approximation is observationally well-motivated due to the azimuthal symmetry (M = 0) in the WMAP-7 
BipoSH measurements in the ecliptic coordinates. All analytical results are verified using corresponding numerical 
simulations on SI maps convolved with NC-beam in the PT-approximation. The numerical simulations also provide 
estimates of the error-bars, hence, a measure of the significance of carious features in the predicted BipoSH spectra. 

As a test case we first work with Elliptical-Gaussian beam with adjustable eccentricity but FWHM that match the 
V and W band WMAP beams. Elliptical-Gaussian beams have clearly interpreted NC parameters and also permits 
a well-controlled perturbation treatment for non-circular deviations (bi m /bio dies rapidly with |m| [1]). This exercise 
also provides a reliable test case to cross-check analytic results with that obtained from numerical simulations of SI 
maps convolved with NC-beams. At the leading order (b^bio), we recover significant CMB BipoSH coefficients Af® 
and Af® 2l with broadly correct qualitative features - relative sign difference in the two spectra and a well understood 
peak structure defined by the acoustic features in underlying C; modulated by spectral shape of 620^/2 • We choose 
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the eccentricity parameters at fixed FWHM to match the amplitude. We conclude that the multiple visits of the 
beam reduces the eccentricity by a factor comparable to that indicated by comparing non-circular beam corrections 
to WMAP Ci from semi-analytic and numerical estimates [2, 40, 42]. 

As may be expected, the Elliptical-Gaussian approximation to WMAP beams, completely fails to recover certain 
qualitative features of the BipoSH measurements, in particular, the change of sign in the V-band Afj 3 at large I. This is 
addressed next by computing the beam-SH for single (A) side WMAP beam maps for V and W bands. The beam maps 
show a marked annular region well beyond a few FWHM from the center where the beam response is negative. This 
leads to the b\2 going negative at large multipoles beyond the beam- width and explains the change in sign of V-band 
BipoSH at larger I. Note that Elliptical-Gaussian can never recover this feature, since bi2 > 0. Again, as expected, 
the amplitude of the CMB-BipoSH based on beam-SH, bi2, from raw beam maps are much higher. An analytic 
treatment of superposition of raw beam transforms to mimic the effect of multiple hits with varying orientation at 
a pixel, leads us to expect a scaling of the leading order term = jibn- However, fi depends on the detailed and 
accurate scan description of each beam and many other finer details of the instruments. Hence we choose to adopt 
a 'phenomenological' approach to determine the effective beam by fitting parametrized /; to the WMAP-7 BipoSH 
spectra measurements assuming it to be entirely sourced by NC-beam effect. A constant fi~am 0.4 provides only 
a fair fit and does not account for the entire signal. A linear correction to the constant scaling, fi = a + 1(3 allows for 
very good fits to the WMAP BipoSH spectra measurements as shown in Fig. 7. 

Complexities of two beam differencing and map- making procedure, deviations from PT-scan approximation, etc., 
may account for the phenomenologically obtained // scaling of raw beam-SH, 6/2- If confirmed in our ongoing numerical 
simulations with realistic WMAP scan and full beam map convolution, we can firmly conclude that NC-beam effects 
arc completely responsible for the WMAP-7 BipoSH measurements. In any case, our results indicate clearly that a 
major portion of the BipoSH signal measured in WMAP-7 may be attributed to an experiment specific systematic effect 
arising from uncorrected NC-Bcam effect. However, it is still possible that the BipoSH measurements are not entirely 
explained by the NC beam effect leaving room for an achromatic component from intrinsic cosmological SI violation. 
Since the WMAP-7 BipoSH signal are strong and significant, other full sky CMB experiments with comparable, or 
better, sensitivity and angular resolution should readily confirm, or disprove this possibility. In near future, good 
quality full sky CMB polarization maps when available can also be studied in the BipoSH representation [32] and will 
provide an independent window in to SI violation phenomena. 
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Appendix A: Beam BipoSH 

Beam-BipoSH are expansion coefficients of the beam response function in BipoSH basis (see Sec III). The most 
general beam-BipoSH in any coordinate system is given by, 



B hh = E ^m^E^'C*) r [ " D l ^ m , (4>,6,p(0,<j>))Yf imi {6,4>) sin 9d9d<f>. 



(Al) 



Wigner-D functions can be expressed in terms of Wigner-rf through following relation, 

D\ nm , (6 9, p) = e"^ d l mm , (9) e- im '" , (A2) 
and reduces to spherical harmonics for m' = 0, 



D l m0 9, p) = y/4n/(2l + l)Y; m (9,4>). (A3) 
In the parallel-transport (PT) scan approximation, the beam orientation, with respect to the local Cartesian 
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coordinate aligned with the spherical (9, 4>) coordinates, does not vary on sky (i.e., p(8, <f)) = p), the beam-BipoSH 

r 2ir 



B hh = V^rY, b ^' e *P~ im ' P E [ exp-C™+"^# f <£ aro ,(l9)< o (0)d(coslSI) (A4) 

= 27r«5MoE & w(^)exp- im '^(-l)™ 2 ^°_ m2i2m2 /W> (A5) 



is derived to be non-zero only for M = 0, using the orthogonality relation, 

r2w 



exp-^ roi+m2 ^# = 27r5 mi ,_ m2 . (A6) 

=o 



We define the notation 



Here we have used d l mm ,(9) = (— l) m ~ m d_ m _ m ,(6). To simplify the analytic expressions, we retain only the leading 
order NC beam spherical harmonic mode m' — 2, assuming mild NC-bcam with discrete even- fold azimuthal symmetry 
where no odd m! modes will contribute. Hence, the summation over m! has three terms corresponding to m! — 0,±2. 
The beam-BipoSH can be then be written as 

ee Btjr + (A8) 

B^ iNC) = 2^0 E ^-™ 2i2 ™ 2 (^-2(5)exp^^_ 2 + 6 i22 (z)exp-^/^ 2 ). (A9) 

First term in Eq.(A8) is the trivial beam-BipoSH, i? ; °; 0<C) , corresponding to the circular symmetric component of the 
beam response function. Non circular part of the beam function ml = ±2, gives rise to beam BipoSH having L =/= 0. 
First, we evaluate the beam-BipoSH due to circular part of beam function. Orthogonality of Wigner-c? functions, 



Ae) = 2TTT' 5 "' (A10) 



implies 



'2 

Using the following property of Clebsch-Gordan , 



hi, -( W 2 \ J(2h + 1) 



(All) 



x;(-i) , ~ m cgSt- m = vw+T)s L0 . (Ai2) 

in 

we obtain 

B^ {c) = V^(-l) h b l20 (z)S lll2 6 L oS M o- (A13) 



Since, 



b l0 (z) = ] /^±^-B l , (A14) 
where Bi is the usual beam transfer function of the circular-symmetrized beam profile, 

Kh (C) = (-1)V2^2 + 1B; 2 <Wlo<Wo. (A15) 
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Next, we analytically calculate the beam-BipoSH P LM< C> due to the non circular part of the beam function, 



B 



lm< nc > 



'In 



(2/i + 1) ^ ,\m ar iL0 

^ dMO^l -1 -' °Zi-m 2 i 2 m 2 



&, 2 _ 2 (z)exp i2 ' / 4 a _ 2 (eX, Wd(cos(9) 

>=0 



+ & (22 (z)exp-^ / dt 22 (0)d l ^ 2O (9)d(cos0) 
Je=o 



(A16) 



In the above expression, the summation is over to 2 . It is convenient to separate the calculation of the to 2 = and 
rest of the to 2 ^ terms. 

Consider the integral for to 2 = 0. For this we use, d l mm , = (— i) m + m d l m , m and expansion of Wigner-rf's in terms of 
associated Legendre polynomials, d l m0 (9) = (—l) m \/{l — m)\/{l + m)!Pr TO (cos 0), and obtain 



=o 



d l %(e)d l j (6)d(cose) 



'(12-2)] r o2 



Pt (cos 9)P h (cos 0)d(cos t 



Using standard recurrence relations of Associated Legendre functions, 

9 ens 

Pf(cos0) = -—-PKcosO) - 1(1 + l)P(cos0), 
sm 

Pi (cos 0) = sin 6P{ (cos 6), 



and orthogonality relations, 

P h (cos 0)P Zl (cos 0)d(cos 9) 



9=0 



2l 2 + 1 



=o 



cos 0P/ 2 (cos 0)P,° (cos 0)d(cos ( 
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if 


(h 


+ h = 





if 


(h 


>h) 





if 


(h 


< h) 


2i 2 + l 


if 


(h 


= h). 



Therefore, for m 2 = 0, the integral simplifies to 



J 0,±2 



-iy 



(2ii + l) 



4tt 



A/ (Z 2 +2)! 



(h-2)\ r 4; 2 2i 2 (; 2 +i) i if 



(i 2 +2)! [(2/ 2 + l) (2Z 2 + 1) 



Next, we consider the m 2 7^ terms in the summation in Eq.(A16). Here, d l ^ 22 
d^ aO (0) using the following recurrence relation, 



sm 



where, 



K 



K. 2 



i| - m| (Z 2 + l) 2 - to 2 . 



(A17) 

(A18) 
(A19) 

(A20) 
(A21) 



(A22) 



h + l 2 = odd) 
h > h) 

h<h) 

h = h). 

can be recursively reduced to 



<fc: o 2 (0) 



K -2«„ 2 



(A23) 



l 2 2 (h + I) 2 ^(4^-1) (Z 2 + 1) 2 (2Z 2 + 1)(2Z 2 + 3)' 
2m 2 



V(z 2 + i) 2 -™ 



+ l)(b + 2)(2Z 2 + 1)' 

v/[a 2 + i) 2 -TO 2 ][g 2 + 2) 2 -TO 2 ] 

+ 1)(Z 2 + 2)(2Z 2 + 1)(2Z 2 + 3) 



= - 2m V!-i)(2z 2 + iy 



and K-2 = 



_ V(^-^)[(z 2 -i)= 



h{h - 1)(4/ 2 - 1) 
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Under reflection Wigner-d's, d* ,(vr - (9) = (-l)'+ m 'd' 



sin 2 6 



Za-l/ 



jZ 2 +2/ 



j/2-2/ 



and then as shown in [2], the integrals for m2 ^ terms are obtained as 



(A24) 
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In general, NC-beams would generate both even-parity and odd-parity beam-BipoSH coefficients 



B 



LM< TO > 



B 



LM<+> 
hh 
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hl 2 



The even-parity beam-BipoSH, 
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if (h = h and l 2 > 2) 



if (h > h and h > 2) 



if (h < h and h > 2) 
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and odd parity-beam-BipoSH, 

' if (h = l 2 ) 

5 M0 [b l22 {z) cxp(-i2p) - b* h2 {z) CX p(*2 /9 )]2^ v /^±li(E| m2 | > o C£?w ima X 

r,,(-5 [f^/ ^t^Mwl 1 :!" 2 ^ + ^ /i^44 '»f ;^- j ^ j ); i) if (z x > z 2 and z 2 > 2) 

r/LM 1 - ' J i\m 2 \ V (Z 2 + l- | m 2 I ) ! ('l + | m 2 I ) ! |m 2 | V (h — 1 — \m 2 \ )!(Zi + |m 2 1 ) U ' v 1 z ^ — I 

^o[^ 2(l)exp(-^)-br 22 (z) cxp( t 2p)]2^ v /^±^E| m 2 | >0 <T_ m2 , im2 x 

r ggi / (i 2 +l-|m 2 |)!(h + l"t 2 ^!" i KK-i / (i 2 -l-|m 2 |)!(ii + |m 2 ^" l N -r/i . . > 9 n 

L|m 2 | V (i 2 + l+|" l2 |)!(i 1 -|m 2 |)! + |m 2 | V -l+|m 2 |)!(Zi -\rn 2 |)! J i 11 W ^ l 2 dllU f 2 ^ Z J 



To avoid any confusion, we reiterate that the above results hold for PT-scan approximation and a NC-beam function 
with discrete even-fold azimuthal symmetry. Other residual symmetries in NC-beam can reduce the set of non-zero 
beam BipoSH further. In particular, if the experimental beam has reflection symmetry, then odd parity beam BipoSH 
will vanish and only even parity ones will be present. This implies that odd parity beam BipoSH can be used as a 
measure of breakdown of reflection symmetry in NC-bcams. 



Appendix B: CMB BipoSH due to Non-circular beam effect 



The cosmological signal in the observed temperature fluctuations is convolved with instrumental beam response 
function. So even if the underlying cosmological temperature fluctuations are statistically isotropic, non circularity of 
the beam can give rise to detections in BipoSH coefficients. In this appendix we provide a detailed derivation of the 
BipoSH coefficient arising from non circular beam presented in Sec. IV. The measured temperature fluctuation map 
AT (hi) is a convolution 

AT (hi) = J dQ n2 B(hi,h 2 )AT(h 2 ). (Bl) 

of the cosmological signal AT(h 2 ) with the beam response function B(hi,n2) that encodes the sensitivity of the 
instrument around the pointing direction, hi and can be expanded in the spherical harmonic (SH) basis, 

AT(h) = Y,aimYi m (h). (B2) 

lm 

Similarly, the cosmological signal decomposed in the SH basis, as 

AT(h)=J2 a imYi m (h). (B3) 
lm 

Beam response function can be expanded in the BipoSH basis, 

B(hi,n 2 )= B hh E Y hmi (fii)Yi imi (ti2). (B4) 

hl 2 LM m\m 2 

Using orthogonality of spherical harmonics, 

J dn A Y lm (h)Y llm ,(h) = (-l) m 's w 5 mm , , (B5) 

we obtain 

AT(ni)=E E (-ir^Ci-m^™!^)- (B6) 

l\m 2 IraLM 

Using the above expansion together with Eq.(B2), we obtain, 

= E (^^^Bhl' 1 Chinl-m- ( B7 ) 
ImLM 
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and the harmonic space covariance as 

/- - \ _ I -l\ m + m ' r, \t>LM TtLM(~iLM riL'M' 

\0-l 1 mi a hni2l — 2^1 ' \ a lm<ll'm') a l 1 l a l 2 V ^m^-m W 2 m 2 2'- 

ImLM I'm'L'M' 

Assuming the cosmological signal to be statistically isotropic, 

{ai m ai' m ') = (—l)CiSu'S m - m i. 
and substituting in Eq.(B8), we obtain the SH-space covariance of the observed map as 

In n \ — ST- (_i\m r i n>LM r>L'M' r LM r L' ' M' 

\ai 1 m 1 ai 2 m 2 1 — \ ' 1 W D hl °Zimi!-m J 2 m2!m 



(B8) 



(B9) 



(BIO) 



ImLM L'W 

As given in Eq.(6), CMB BipoSH coefficients are related to the SH space covariance matrix in Eq. (B8), leading to 

E/-1 r>LM tjL'M' \ y ( -i\m/~iLM /~iL' M 1 s~<LiMi fail} 
^l D hl n l 2 l Z^i \ l ) u lim 1 l-m L/ hm 2 lm L/ l 1 mihm2 \ all J 



hh 



ILL' MM' mm 1 m 2 

The sum over product of three Clebsch-Gordan coefficients can be written compactly in terms of a 6-j symbol, as 



£(-l)-"C^C3^C# _ a = A'^C, 

apS cf 



a b c 
c tJV1 e / d( ' 



(B12) 



where Ky = (-l) h + c + rf +/ and ]J cf = ^{2c + 1)(2/ + 1). In a PT-scan approximation, M = 0,M' = 0, Mi = 0. 
Hence, we obtain the expression in Eq. (32) for CMB BipoSH coefficient from NC-beam with PT-scan approximation 

A hh h = <W £ C l BftBff(-l)W- L * ^(2L + l)(2L> + l)C L L &J! £, f\ . (B13) 

ILL' { J 

The Clebsch-Gordan coefficient C^Q° L , a is zero when the sum L + V + L\ is odd valued, hence, enforces the condition 
that the summation in the above expression is limited to L + V + L\ being even-valued. When the beam function 
has an even fold azimuthal symmtery and reflection symmetry beam-BipoSH coefficients are restricted to even parity 
and follows l\ + 12 = even, then L and L' are restricted to even multipole values. Thereafter, due to the presence of 
Clol'O' -^i takes up even multipole values. 



Appendix C: Effective beam averaged over multiple hits with varying orientations 

We present an analytic treatment to cover NC-beam effect incorporating the multiple hits at any pixel n by the 
NC-beam with fixed shape but at varying orientations, Pj(n). This information can be obtained from the instrument 
design description and the scan-strategy of the experiment over the duration of the data acquisition. 

The observed temperature anisotropy for a single hit by the beam at a direction n = (6, <f>) with orientation, Pj(n) 
is given by 



AT(n) = / dn n ,B{n, ti- Pl {h)) AT(n') , (CI) 

where AT(h') is the cosmological signal. Simply generalizing the notations for Eqs. (9) to (11) in Sec. Ill, we write 
B{n,h'; Pi {h)) = 2~2b bn (h; Pl )Y lm (h r ) (C2) 



= E 

lm 

= E 



lm L m 



lm L m 



2^2b lml (z)D l mm ,((i),e,p t (h)) 

. m' 



Y hn (h') 
e- im *Y lm (n') 



(C3) 
(CM) 
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where we relate the beam function in terms of beam-SH, bi m i(z) as done in Eq. (11). 

Next, we define an effective beam-SH 6^,(z) averaged over N(h) multiple hits at pixel, ft, and over the whole sky 



J .7 = 1 ^ > m> 



bi m '(z) (C5) 



where it is easy to show that 



I 111 



N(n) 

= ]T J dn n d l mm , (6) e-™* e" w « W . (C6) 



The above provides the general expression for the beam-SH of the effective beam for a single beam experiment given 
the beam-SH from the raw beam map and entire scan pattern history, pj(n) . Note, it is tedious, but conceptually 
straightforward, to extend this formalism to account for two beam differencing even with different A and B side NC- 
beam by separating the averages over b^*(z; pj ' B (ft) with appropriate signs and weights. However, the limitation 
of this approach is that it does not capture the effect of the linear algebraic step of map-making carried out to covert 
difference measurements between different pixels into temperature fluctuations map. 

In this paper, we have evaluated analytic expressions only up to to the leading order effect at m = 2. Since the 
WMAP beams seem mildly NC-beams where \bi m \/bio become smaller with larger \m\ (at least for lowest few values 
of m > 2 that were explicitly checked). Hence, we provide explicit expressions for bf 2 from Eq (C6), limiting to 
beam-SH bi m {z) = 6/o<5 m o + bu5 m i + bi 2 S m2 (i- e -> assuming mildly NC-beam ) as 

b^(z) = bioWFio + bnWFn + buWFu (C7) 

Fio 



J dn n e l2(t, d l 20 = J dQ n Y 2Q (n) = b 00 5 l0 (C8) 

Fn = I dO n e- i2 <* [4-i(% iKft) +4i(0)e" ip(A) ( C9 ) 
"* i=l 

r N{A) 

Fn = / dfi„c- i2 * J2 [d l 2-2^y 2p{A) + d l 22 (e)c- l2p(A) . (CIO) 



The first term integrates to zero for I > 0. In our analytic expressions, we ignored, 6;i(ft) , since it drops out in the 
PT-scan approximation even without reflection symmetry assumption. This may not necessarily be non-negligible 
here due to deviations from PT-scan approximation but we ignore it for simplicity. We finally obtain the scaling used 
in Sec. VB 

b$(z) = F l2 b l2 (z) 

N(n 



F i2 = J dnne-^ ^ [4 _ a (0) c l2p ^' l) + d l 22 {9) c^ 2 "^] . (Cll) 



Invoking the approximation that Pj(h) « pj is constant for each pass of the beam on the sky and equal number 
of visits, fi = a, can shown to be a constant, \a\ < 1. We use a simple model of the WMAP scan strategy [44] that 
provides an estimate of a ~ 0.4 for identical A & B side beams at a relative angle of ~ 140° consistent with the 
best-fit a obtained from the BipoSH spectra. The simple WMAP scan model ignores, flagged portions of the actual 
time-line data, finer details of the location of the different frequency feed-horns from the spin-axis on the focal plane, 
etc., that may, or may not, be important. In the absence of unambiguous means to determine fi, we resort to making 
a polynomial (in I) fit to /; and find that a linear term is sufficient to provide good fits to the WMAP-7 BipoSH 
spectra measurements. From the analytic expression in Eq. (Cll), apprpriate non-polynomial parametrized forms of 
fi may yet be found to be better motivated. 
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Odd parity BipoSH coefficients are a measure of the breakdown of reflection symmetry in instrumental beams. Expressions 
for the odd parity BipoSH coefficients generated in absence of reflection-symmetry of the beam are provided in the appendix 
for completeness. 

[46] Departure from reflection symmetry in the beam of a full-sky CMB experiment, if ignored, also causes leakage of power from 
the ~ 500 times stronger CMB dipole signal into higher multipole, most importantly, contaminating the CMB quadrupole 
moment of the angular power spectrum. This has been studied and estimates on WMAP beam maps indicates the effect 
of reflection breakdown symmetry is expected to be small, but not negligible [44]. 

[47] Note that this factor in WMAP-BipoSH estimator strictly restricts BipoSH considerations to the even parity sector since 
Cio?'o = f° r °dd values of the sum L = l+l' . In the context of NC-beams, this would be a handicap if reflection symmetry is 
violated leading to odd-parity BipoSH coefficients (see Sec. III). Also it is blind to a number of other interesting possibilities 
with odd-BipoSH signals. 



